Current fluctuations near to the 2D superconductor-insulator quantum critical point 
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Systems near to quantum critical points show universal scaling in their response functions. We 
consider whether this scaling is reflected in their fluctuations; namely in current-noise. Naive scaling 
predicts low-temperature Johnson noise crossing over to noise power oc £ z /( z+1 ) at strong electric 
fields. We study this crossover in the metallic state at the 2d z=l superconductor/insulator quantum 
critical point. Using a Boltzmann-Langevin approach within a f /TV-expansion, we show that the 
current noise obeys a scaling form Sj = T$[T/T e ff(E)] with T a ff oc y/E. We recover Johnson noise 
in thermal equilibrium and Sj oc \[E at strong electric fields. The suppression from free carrier shot 
noise is due to strong correlations at the critical point. We discuss its interpretation in terms of a 
diverging carrier charge oc 1/V~E or as out-of-equilibrium Johnson noise with effective temperature 
oc y/E. 

PACS numbers: 



The notion of quantum criticality provides one of the 
few general approaches to the study of strongly correlated 
quantum many-body systems 1]. The scale invariance 
that characterizes the zero temperature critical point 
leads to characteristic, universal power-law dependences 
for various quantities in its proximity; these dependences 
can be computed within a continuum field theory. A 
number of theoretical works have calculated the impact 
of such universality upon conductivity 0,0,0. These cal- 
culations provide robust experimental predictions. Two 
recent works of Dalidovch and Philips[5j, and Green and 
Sondhi 0], have extended the analysis to see whether uni- 
versality persists out of equilibrium. Surprisingly it does; 
at least in transport. 

In this work we consider another experimentally mea- 
surable quantity: current fluctuations. Measurements 
of current fluctuations are usually restricted to meso- 
scopic samples in order that the relative fluctuations of 
current be significant. Mesoscopic samples have rela- 
tively few conducting channels giving a total conduc- 
tance of order the conductance quantum e 2 jh. Quantum- 
critical systems may also satisfy this criterion. For ex- 
ample, at the two-dimensional superconductor-insulator 
(Sl)transition the conductance is also of order the con- 
ductance quantum^, 0, HI El It would be revealing, 
therefore, to compare the current fluctuations at the 
metallic critical point of bosons with the more familiar re- 
sults for non-equilibrium current noise in diffusive Fermi 
liquids 0. Transport properties are readily measurable 
and theoretical predictions regarding their scaling and 
universality at two-dimensional quantum critical points 
(in charged systems such as SI, and Quantum Hall Effect 
transitions) have stimulated many experiments. Non- 
equilibrium current fluctuations may provide an indepen- 
dent window upon these quantum critical points. 

There are two contributions to current noise; John- 



son noise and Shot noise. Since it is a reflection of 
the fluctuation-dissipation relation, Johnson noise is ex- 
pected to be unaffected by correlations between carriers. 
On the other hand, shot noise reflects the fact that charge 
is carried in quanta. If there are strong correlations be- 
tween carriers, we expect it to be significantly altered. 

The situation in a non-equilibrium steady state is trick- 
ier; Johnson and shot noise cannot easily be separated. 
On the one hand, one may identify an effective temper- 
ature scale T e ff for the out-of-equilibrium distribution 
and treat the current fluctuations as Johnson-like. Di- 
mensional analysis suggests T e ff oc E z ^ 1+Z ^ and current 
noise proportional to E z ^ 1+Z \ On the other hand, one 
may identify the out-of-equilibrium noise with shot noise 
and a carrier charge that diverges as E~ 1 ^ 1+z \ 

We present a microscopic calculation that recovers 
these results for the superfluid to insulator transition of 
two-dimensional bosons with particle-hole symmetry. We 
use a Boltzmann-Langevin |f 1 1 1 ( I] approach within a 1/N- 
expansionjj] to analyse current noise near this transition. 
This "kinetic theory of fluctuations" approach is known 
to give the correct result for the leading noise correla- 
tions in the diffusive metal, but needs to be modified 
for higher moments jjjj. The system is considered to 
be in thermal contact with a substrate. This, in con- 
junction with the internal scattering in the system, leads 
to a uniform, non-equilibrium, steady-state distribution, 
provided that the length of the sample is greater than 
the correlation length 6] . This situation is quite different 
from that considered elsewhere, where the only coupling 
to the outside world is through the leads. It permits the 
current-noise properties to be entirely determined by the 
internal scattering and so to be universal. 

Our main results are as follows: in thermal equilib- 
rium, we recover Johnson noise with noise power Sj = 
AcrT (as we must in order to satisfy the equilibrium flue- 
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tuations dissipation relation). This crosses over at very 
large electric fields to Sj oc y/~E. This is strongly sup- 
pressed from the usual shot noise result of Sj oc j oc E 
due to the strong correlations at the quantum critical 
point. This large field result may be considered as either 
shot noise with a carrier charge which diverges as 1/ \f~E~ 
or alternatively as a non-equilibrium Johnson noise with 
effective temperature T e ff — \fficeEir 2 / Ak. Between 
these two limits, we expect the current noise to obey 
a universal scaling form Sj = T$[T e ff(E)/T]. 

We begin by outlining the field theoretical formulation 
of the symmetric, bosonic, superfluid to insulator transi- 
tion. We then give Boltzmann and Boltzmann-Langevin 
equations for the system in thermal equilibrium. The 
latter is used to deduce the Johnson noise. Next, we 
turn to the Boltzmann and Boltzmann-Langevin equa- 
tions for the zero-temperature system under a finite elec- 
tric field. This is used to deduce the out-of-equilibrium 
current-fluctuations. Finally we turn to a discussion of 
the implications of these results. 

Field Theory The critical region of the symmetric 
superfluid to Mott insulator transition phase diagram 
is described by a charged scalar field with a quartic 
interaction[lj,|2|: 
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J d d x [irn + v</> f V0 + mW + ao^) 2 ] , (i) 



where <j) is the complex scalar field and n is its conju- 
gate momentum. These satisfy the usual commutation 
relations [</>(x, t), n(y, t)] = i<5(x — y). It is convenient to 
choose the bare interaction A to have its fixed point value 
u* A 3_d with momentum cutoff A, pj although we will not 
need the precise, regularization-dependent, value. At the 
zero-temperature critical point, the renormalized mass is 
zero which corresponds to a particular choice m* of the 
bare mass. The effects of applying an electric field, E, 
are included by minimally coupling to a vector poten- 
tial; Vcj) — > D<j> = (V + ieA/%) (j). We choose the gauge 
A = Et. This is equivalent to a contact-free measure- 
ment (for example by placing the system on a cylinder 
and inducing an EMF on its surface by uniformly in- 
creasing the flux through it) . We believe that our results 
carry over unchanged to the case with good contacts at 
the ends. 

The normal modes of this Hamiltonian (in the absence 
of interaction) are charge density fluctuations. These oc- 
cur with positive and negative charges, corresponding 
to decrease or increase in charge density from the av- 
erage. We determine the current fluctuations near to 
the critical point by considering a Boltzmann equation 
for the occupation of these modes. We use a' ftf and 
a/b to represent the creation and annihilation of posi- 
tively/negatively charged density fluctuations. 

The Boltzmann equation in thermal equilibrium and 
within a l/A^-expansion of the Hamiltonian Eq.Q is 



given byH 

[8t + [eE/h).dk] h 

7kq/k (1 + /q) - 7qk/q (1 + /k) 
7kq/k/-q ~ 7qk (1 + /-q) (1 + /k) 

7k, q and 7k, q describe the matrix elements for 
particle-particle and particle-anti-particle scattering, 
respectively ^2 an d 



= - dq 



•(2) 



/(k,i) 



dq(al 



q/2 



(*)o k _q/ 2 (t)) 



denotes the distribution function, with a similar distribu- 
tion /(k, t) for the negatively charged modes. Particle- 
hole symmetry implies the relation /(k, t) — /(— k, t). 
We restrict our explicit consideration to the positively 
charged channel; final expressions include the negatively 
charged channel through appropriate factors of two. 

The Boltzmann-Langevin equation is an equation de- 
scribing the stochastic evolution of fluctuations, <5/ q (r, t), 
in the distribution function about its equilibrium or 
steady-state. In the present case, it is given by 

[d t + v.d r + (eE/7i).d q ] Sf k = 7k [Sf k - <5/ k ] + r? q , 



dq 
dq 



7k q i | £ o/k 



7kq 



l + /k 

/-q 
l + /k 



7qk 1 ; , 0/q 



5fk~ Tqk 



l+/q 
l+/k 
/-q 



7?q 



V-q 



(3) 



where 



M kq 

7k 



dqM kq (5/ q 



7qk 1 + /k 7-gk 1 + /k 

7k 1 + /q 7k /q 

/ rfq ( 7kq iTA + ^T 



+ h 



(4) 



The position and time labels of <5/k(r, t) and ?7 q (r, t) have 
been suppressed for compactness and the detailed bal- 
ance conditions 7q,k/q(l + /k) = 7k, q /k(l + /q) and 
7 q ,k(l + /- q )(l + /k) = 7k,q/k/- q have been used in 
order to simplify the right-hand-side. The time and mo- 
mentum derivatives on the left-hand-side of Eq. @ follow 
directly from the Boltzmann equation @ ■ The first term 
on the right hand side is the linearised form of the scat- 
tering integral. The second is a stochastic term describ- 
ing fluctuations in occupation number. It is determined 
by assuming that the scattering processes are indepen- 
dently Poisson distributed. Under this assumption, the 
quadratic correlations of 7y q (r,t) are given by 

(T^M^rV)) 



nptu 

m 



= 2S rr 'St 



(27T) 2 £ q , q , / dk 7qk / q (l + /k) 
-7qq'/q(! + /q') 

+ (27r) 2 <5 qiq , J dk7 qk / q /-k 

— 7qq'/q/— q' 
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2tfr,r'$t,t'7q'/q(l + /q)[<W ~ M q'q] 



(5) 



where the detailed balance conditions and Eq. Q have 
been used in order to simplify these expressions. The 
factor of two comes from the contributions of in- and 
out-scattering processes. 

Current fluctuations may be calculated using the 
Boltzmann-Langevin equation as follows: first we deter- 
mine the local fluctuations in occupation. In the limit 
of long times and long length-scales, we may neglect the 
derivative terms on the left-hand-side of Eq. 0. The 
fluctuation in the distribution function is then given by 

<5/k - Sfu = Vk/lk- (6) 
This integral equation has a formal solution given by 

*/ k (x, t)= J dq[l - M^r^x, t)/ 7q , (7) 

where we are using a matrix notation for functions of 
momentum with [l]kq = (27r) 2 <5(k — q) and the inverse 
of a matrix N denned as / dq[N]^ [N] q k' = [l]kk'- 
The resulting fluctuation in current is given by: 



Sj= dpvpSfq, 



(8) 



where Vk = <9k£k- Using the noise correlations from 
Eq. JSJ), after some algebra, we obtain the following ex- 
pression for the correlation of current fluctuations in ther- 
mal equilibrium |l4j |: 

(6j a (r,t)S (r',t')) = 5 a ,05 TtT iSt,t> J dpdqv p .v q 

x[l-M] pq / q (l + / q )/ 7q (9) 

This result should be compared with the result for 
fermions0 where the Bose enhancement factor / + 1 is 
replaced by a Fermi factor 1 — /. The low frequency 
current noise Sj can now be expressed as: 



dfcdr<j(r,i)j(0,0)) 



(10) 



In Thermal equilibrium, the current fluctuations re- 
duce to the Johnson result Sj = Aaj-T, with the uni- 
versal conductivity cjt determined by a linear response 
expansion of the Boltzmann equation (J2J) 1(|. As the 
electric field is increased away from zero, symmetry un- 
der E — > — E requires that the lowest order correction to 
this result must be proportional to E 2 . A closed form 
expression for this correction may be obtained by includ- 
ing the lowest order, static distortion of the distribution 
function by the electric field, Sf s Hi . The resulting ex- 
pression is cumbersome and not very revealing. We note, 
however, that the corrections to the distribution func- 
tion are proportional to powers of TiceE/{kT) 2 (to see 
this, rescale the momentum integral in the expression for 



6f s in and use the fact that for critical, relativistic 
bosons |vk| = c). The current noise is, therefore, ex- 
pected to be a universal scaling function of HceE/ (kT) 2 - 



Sj = Tf [hceE/(kT) 2 ] = T<f>[T eff (E)/T], 



(11) 



where we have defined T eff (E) = VhceETT 2 /4k (the rea- 
son for the numerical factor in this expression will be- 
come apparent later). The low-temperature limit of this 
scaling function is $(0) = 4<tt, where o~t — Q-21^Ae 2 /h 
given by the universal conductance found in thermal 
equilibrium^, 0]. Next, we turn to a calculation of the 
high-field limit of this scaling function. 

Out-of-Equilibrium. Under the application of a strong 
electric field, the system is driven far from thermal equi- 
librium. In this situation, the Boltzmann equation 
misses important field-induced tunneling processes (pre- 
cisely analogous to Zener breakdown or the Schwinger 
mechanism). A Boltzmann equation for this situation 
was previously derived within a 1 /N-expansion in Ref. ; 



[d t + (eE/fi).^] /(k) 



J(k) + e -* nce f eE 5 ( 



nk.E 



eE 2 



(12) 



where to lowest order in 1/N|lq. e k = k 2 an d 



rk(t) = 8c^m/^vi k 'i+*. 



N e k 
2 1 



£k' 



/(k',t) 



(eE/hf^c 1 '*. 



(13) 



Eci. (|12|) may be integrated to obtain an explicit solution 
for the out-of equilibrium distribution function^, IT^ | . 

Following the same procedure as previously, we may 
derive a Boltzmann-Langevin equation describing fluctu- 
ations in occupation number. This has a simpler form 
than in thermal equilibrium since there is just one pro- 
cess that dominates the contribution to the fluctuations 
in occupation number in this strongly out-of-equilibrium 
limit. The result is (suppressing position and time labels) 

[d t + v.a r + E.a k ] sf k = -r k <5/ k + % 

(?7 q (r, t)rks (v',t')) = (2ir) 2 5 ry 5 t , t ,5(q - q')r q / q (14) 

Solving, as before, in the limit of long times and large 
distances, the fluctuation in occupation number is given 
by 5fq = r/q/Tq and the correlation of fluctuations in 
current are given by 

(j a {r,t)j (r',t')) 

= 2e 2 5 a>/3 f dpdq^^<»fc(r > f)» ?q (r' > f')> 

J 1 pi q 

= 2e 2 c 2 5 a , f) 6 r y5 t , t , J dp^, (15) 
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where we have used the fact that v p = c 2 for relativistic 
bosons in the large N limit and the factor of two comes 
from the contribution of particles and holes. Using the 
explicit form of the non-equilibrium distribution function 
and performing the integrals, one obtains the high-field 
contribution to current noise as 



Ua(r,t)j/3(r ,t)) = 5 a> p5 t ,v6 T ,r> I — 



ecE 



16) 



The large field current noise (|16[) is dramatically re- 
duced compared with the linear E-field dependence of 
shot noise expected for uncorrelated charge carriers. It 
is tempting to ascribe this to an effective charge of carri- 
ers near to the quantum critical point proportional to 
1/VE. This effective charge shows a dramatic diver- 
gence as one approaches the quantum critical point re- 
flecting the strong correlations of the system. Whether 
this picture can be carried through requires an analysis 
of higher order statistics of the current noise^l|. An a l~ 
ternative, and perhaps more revealing, interpretation of 
Eq. (|lbl is as a non-equilibrium equivalent of Johnson 
noise. The non-equilibrium analogue of the fluctuation 
dissipation relation is recovered if we identify an effective 
temperature T e ff = \JhceE it 2 /Ak. This is deduced from 
4a E T ef f = 2{Nir/8) 2 y/ecE/ft ^ where a E = {Nn/8)e 2 /h 
is the large-field conductivity [6j. The current-noise scal- 
ing function has high-field asymptote $[x — » oo] = Ao~ex. 
To estimate the magnitude of this effect in a physical set- 
ting such as the SI transition in MoGe thin films 0] , we 
can take E — 0.5V/m, and estimate c ~ vp = 10 6 m/s. 
We obtain T e // ~ 400mK, which is an order of magni- 
tude larger than experimentally accessible temperatures. 
The current noise will, therefore, be enhanced over its 
thermal value. 

In conclusion, we have considered the current noise at the 
2-dimensional, z=l superconductor to insulator quan- 
tum critical point. We find that this noise follows a 
universal scaling function, Sj = T$[T e / f(E)/T], with 
T e ff = v 'TiceE-rr 2 1 '4. This scaling function recovers John- 
son noise in thermal equilibrium and crosses over to 
an unusual -\/i?-dependent non-linear shot noise or non- 
equilibrium Johnson noise at high fields. This is a par- 
ticular case of a more general E z '( 1+Z > scaling expected 
for high-field current noise. In this way, current-noise 
may reveal the universal non-linear scaling exponents 
predicted near to quantum phase transitions. 
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7k,. 



2 

N 



and 



7k,. 



2 

N 



dtt Imn _1 (q- k, fi) 
ty 4e k e q 

" [1 + n(fi)]£(ek - e q - Q) 
+n(£i))5(e k - e q + O) 



dfiTmn _1 (k + q,n) 
7T 4e k e q 
x [1 + n(n)]<5 (e k + e q - Q) 



Il(q, SI) is the polarization loop for particles or anti- 
particles and is determined by the critical scattering; we 
will give no further details of it here except to use its crit- 
ical form in our final calculations. For a complete com- 
parison with the expressions given in Ref.Q, one must 
make use of the identity n(fi) = e~ n/T (l + n(fi)). 
[13] In fact, there is a subtlty here. Particle-hole scattering in- 
duces correlations between noise in the particle and hole 
channels. Eq. JSJ incorporates these effects such that we 
may treat particle and hole scattering as if they were in- 
dependent obtaining the correct result. Strictly, the noise 
correlations are given by (suppressing space and time la- 
bels) 



(2vr) 2 5 q , q , J dk 7q , k / q (l + / k ) 
-7q,q'/q(l + /q') 

+ (27r) 2 <5 q>q , J dk7 qik /i, /g /_ k _ 

(^V) = 2<5 ri ^5 M / [7q,q'/q/-q' ] 



writing the total contribution to the current fluctuation 

as 



<5j = / rfpv 



7p 

with these noise correlators we recover the results derived 
in the text. 
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[14] In matrix notation the equation reads: 

<77 q (r, ijife, (r', t')) = 2* r , P , 5 4 , t , 7q , / q (l + / q )[l - M] q , q 

Combining this with Eq. (J7J we find the fluctuations of 
the distribution function are correlated as: 



<a/p(r,t)<5/ q (r',i')> = 25 r y S t>t ,[l-M] 



-l/q(l + /q) 



hence correlators of the current fluctuations © are: 

{Sj a (T,t)6f > {i f ,t')) = 8 a ,p6 tit ,8t,t, f [l-M] pq / q (l + / q )/ 7q 

[15] Within the 1/N-expansion for the out-of-equilibrium sys- 
tem, the gap does not contribute to the current response 
at lowest order in 1/N and is neglected here (See also 
i). 

[16] The linear response conductivity may be calculated from 
the Boltzmann equation @ in the usual way. Expanding 
to linear order in the applied electric field, we find 



J 7k 

where 8f now indicates the static, linear distortion of 
the distribution from the thermal distribution / in re- 
sponse to the electric field. We may formally invert this 



integral with the result 

Sfi = J dq[l-M]^^E.ft,/q 

Finally, using j = Jdkv k <5/k, and 9 q / q = v q (l + 
/q)/q/feT we may deduce the following linear response 
solution for the conductivity in thermal equilibrium: 

kTa = 2 / dpdqVp.v q [l-M] pq / q (l + / q )/ 7q 



the factor of 2 arising from the contribution of both par- 
ticles and anti-particles to the conductivity. This may be 
compared against Eq.@. When calculated within a 1/N- 
expansion, the result of carrying out these integrals is a 
universal conductivity in thermal equilibrium, given by 
ot = N0.1077e* 2 /h\$. 

[17] The explicit solution for the distribution function is / k = 
O(K^) exp [— -kK\ — 7«/2K"|| I , where the momenta have 
been rescaled according to K = \t/y/ eE/hc and 7 is 
determined by T k = ■y(eE/h) 3/i c 1/i / v /2fcjj". In terms of 
these variables, Eq. il 1 2li may be integrated to obtain the 
above solution. The integrals in Eq. 1131 may be carried 
out noting that ku 3> k± to obtain r y 2 — 8/Ntv 2 . 

[18] N. Mason and A. Kapitulnik, Phs. Rev. Lett. 82, 5341 
(1999). 



